LTH 516 



O 
O 



c^ 



Is the up-quark massless? 



UKQCD Collaboration 



A.C. Irving, C. McNeile, C. Michael, K.J. Sharkey and H. Wittigj] 

Division of Theoretical Physics, Department of Mathematical Sciences, University of 



^ ' Liverpool, Liverpool L69 3BX, UK 

<' 
(N 

(N 
> 

m 
(N 

O 

o 



Abstract 

We report on determinations of the low-energy constants 05 and ag in the ef- 
Qh. fective chiral Lagrangian at O(p^), using lattice simulations with A'^f = 2 flavours 

r~| , of dynamical quarks. Precise knowledge of these constants is required to test the 

hypothesis whether or not the up-quark is massless. Our results are obtained by 
K/i I studying the quark mass dependence of suitably defined ratios of pseudoscalar me- 

H ' son masses and matrix elements. Although comparisons with an earlier study in 

the quenched approximation reveal small qualitative differences in the quark mass 
behaviour, numerical estimates for a^ and as show only a weak dependence on 
the number of dynamical quark flavours. Our results disfavour the possibility of a 
massless up-quark, provided that the quark mass dependence in the physical three- 
flavour case is not fundamentally different from the two-flavour case studied here. 
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1 Introduction 

A massless up-quark represents a simple and elegant solution to the strong CP problem. 
Consequently, the question of whether or not rriu is indeed zero has been the subject 
of much debate over many years (for a review see ref. |l[|). Traditionally the problem 
is studied in the framework of Chiral Perturbation Theory (ChPT). Although the most 
recent estimates point to a non-zero value for the ratio m^/m^ [Q, the situation is com- 
plicated by the presence of a hidden symmetry in the effective chiral Lagrangian 0] . This 
so-called "Kaplan-Manohar ambiguity" implies that m^/rrid can only be constrained after 
supplementing ChPT with additional theoretical assumptions. Although the validity of 
the commonly used assumptions is plausible @,||, it is clear that the question whether 
the up-quark is massless cannot be studied from first principles in ChPT. 

More recently attention has focussed on lattice simulations to tackle this problem. A 
reliable, direct lattice calculation of m^, however, presents considerable difficulties, even 
on today's massively parallel computers. It was therefore proposed to use a more indirect 
approach, based on a combination of lattice QCD and ChPT P-§]. The aim of this method 
is a lattice determination of the so-called "low-energy constants" in the effective chiral 
Lagrangian, whose precise values are required to constrain mu/wid using chiral sjTiimetry 
and phenomenological input. A variant of this proposal, which allows for a determination 
of the low-energy constants with good statistical accuracy was discussed in ref. [|] and 
tested in the quenched approximation. 

Here we extend the study of ^ to QCD with two flavours of dynamical quarks. While 
this addresses the important issue of dynamical quark effects, it still does not correspond 
to the physical three-flavour case, and thus we are yet unable to give a final answer to 
the question in the title of this paper. Nevertheless, our study represents an important 
step in an ultimately realistic treatment of the problem, by studying the dependence of 
the low-energy constants on the number of fiavours. If our results can be taken over to 
the physical case without large modifications - and there are indications that this is not 
unreasonable - then the possibility of a massless up-quark is strongly disfavoured. 

In Section ^ we briefiy review the Kaplan-Manohar ambiguity and its relevance for the 
value of rriu. Section || contains details of our lattice simulations, whose results are de- 
scribed in Section |^. In Section |^ we discuss the implications of our findings and present 
an outlook to future work. 



2 Low-energy constants and m^ = 

In order to make this paper self-contained, we briefiy review the implications of the 
Kaplan-Manohar ambiguity for the ratio m^/md. The strategy to address the problem in 
lattice simulations will then become clear. A more complete discussion can be found in 



refs. 0,|,0,i. 

A determination of m-a/m^ in ChPT which is able to distinguish between a massless and a 
massive up-quark requires precise knowledge of the first-order mass correction term Am- 
At order p^ in the chiral Lagrangian it is given by [p!0|-p!2[] 

2 2 

Am = ,. p \2 *^2"8 ~ "5) + chiral logs, (1) 

where F^ = 93.3 MeV is the pion decay constant, and 0:5, as are low-energy constants, 
whose values have to be determined from phenomenology. Throughout this paper we 
adopt a convention in which the low-energy constants a, are related to the corresponding 



constants Lj of ref. [11| through a^ = 8(47r)^Lj. Furthermore, we always quote low-energy 
constants in the MS-scheme at scale /x = At^Ft^. 

The value of as can be extracted from the ratio of pseudoscalar decay constants, Fk/Ftt, 
and is given by 

as = 0.5 ±0.6. (2) 

By contrast, there is no direct phenomenological information on ag or the combination 
(2a8 — as). Although as is contained in the correction to the Gell-Mann-Okubo formula, 

i.e. 

2 2 

yoo ^ 777 

Agmo = .^ p .2 ("5 ~" 12a7 - Gas) + chiral logs, (3) 

its determination requires prior knowledge of aj. It is at this point that the Kaplan- 
Manohar (KM) ambiguity becomes important. It arises from the observation that a 
simultaneous transformation of the quark masses 

m^ -^ m^ + Xrridms, rrid -^ rrid + XrUsm^, rUs -^ rng + Xm^rrid, (4) 

and coupling constants according to 

ae^ae + X^-— — , ay -> ay + A— — — , as -^ ag - A^— — , (5) 

4-Do 41^0 ^-dq 

leaves the effective chiral Lagrangian invariant. Here, A is an arbitrary parameter, and Fq, 
Bq are coupling constants in the lowest-order chiral Lagrangian.^ Thus, chiral symmetry 
cannot distinguish between different sets of quark masses and coupling constants, which 
are related through eqs. (^) and (|). Indeed the correction Agmo is invariant under the 
above transformations. The value of ay can be fixed by invoking additional theoretical 
assumptions, such as the validity of large- A'"c arguments for N^. = 3. In accordance with 
these assumptions, Leutwyler [Q constrained the correction term Am to be small and 
positive: 

< Am < 0.13. (6) 

'^Fq coincides with F^r at lowest order. 



This gives 771^/111,^ = 0.553(43) and hence a non-zero value for rriu. The "standard" values 



for a^ and ag which are compatible with eq. (y) are ||T0Hl2| 

07 = -0.5 ±0.25, as = 0.76 ±0.4. (7) 

In view of the importance of the strong CP problem, one may regard any analysis based 
on theoretical assumptions beyond chiral symmetry as insufficient. In particular, since the 
uncertainties in the estimates for ai and as are quite large, the possibility that m^ = 
does not appear to be ruled out completely. A massless up-quark would require 

a7 = 0.25 ±0.25, ag = -0.9 ± 0.4, (8) 

resulting in a large, negative first-order correction Am- In order to decide which scenario 
is realised and to pin down the value of m-a one has to replace the theoretical assumptions 
by a solution of the underlying theory of QCD. 

The KM ambiguity implies that the low-energy constants aj, as (and ag) can be deter- 
mined from chiral symmetry and phenomenology only if the physical quark masses are 
known already. At this point it is important to realise that QCD is not affiicted with the 
KM ambiguity, and that the formalism of ChPT also holds for unphysical quark masses. 
Since quark masses are input parameters in lattice simulations of QCD, their relations to 
hadronic observables need not be known a priori. Hence, the low-energy constants can be 
determined by studying pseudoscalar meson masses and matrix elements for unphysical 
quark masses and fitting their quark mass dependence to the expressions found in ChPT. 
In this way it is possible to determine as and - more importantly - the combination 
(2a8 — as) directly in lattice simulations. 

3 Lattice setup and simulation details 

In ref. ||^ it was shown how the low-energy constants can be extracted from lattice data 
for suitably defined ratios of pseudoscalar masses and matrix elements, Rm and Rp. Here 
we repeat their definitions in order to explain the necessary notation. For more details 
we refer the reader to the original paper 0. 

The actual determination of the low-energy constants proceeds by studying the mass 
dependence of -Rm and Rp around some reference quark mass niy-cf. As pointed out 
before, rrirei does not have to coincide with a physical quark mass |T^, as long as it is 
small enough for ChPT to be applicable. 

In this paper we will be concerned with simulations of "partially quenched" QCD, where 
valence and sea quarks can have different masses. Let us therefore consider a pseudoscalar 
meson with valence quark masses mi and m2 at a fixed value of the sea quark mass m^^^. 
In [PI the dimensionless mass parameters 



Bo(mi + m2) 2BQm 



ref 



y= (4vrFo)^ ' y-' = j4^^^ ^ = y/yreu (9) 



were defined, as well as the ratios 




Ry{x) = Fps{y)/Fps{yrei), x = y/y,,i. (10) 

Here Fps is the pseudoscalar decay constant, and Gps is the matrix element of the pseu- 
doscalar density between a pseudoscalar state and the vacuum. The parameter x denotes 
the fraction of the reference quark mass at which the ratios -Rm and Rp are considered. 

As emphasised in 0, Rm and Rp are well-suited to extract the low-energy constants, 
since ratios are usually obtained with high statistical accuracy in numerical simulations. 
Furthermore, any renormalisation factors associated with Fps and Gps drop out, so that 
Rm and Rp can be readily extrapolated to the continuum for every fixed value of x. As 
has been shown in the quenched approximation ^, discretisation errors in Rm and Rp are 
very small, so that good control over lattice artefacts is achieved. This issue is important, 
since the determination of low-energy constants requires an unambiguous separation of 
the mass dependence from effects of non-zero lattice spacing. 

We now give the expressions for Rm and Rp in partially quenched QCD which are relevant 
for our study. They were obtained using the results of ref. ^^ and are listed in appendix A 
of 1^. Our reference point y^-^f was always defined at 



mi = ma = m"""" = m^f, 2/ref = 77-77^. (H) 



2 Bom j-cf 

(47rFo) 



In order to map out the mass dependence of Rm and Rp for a fixed value of m^'^^, we 
have considered the cases labelled "VV" and "VSl" in ref. [0. The first uses degenerate 
valence quarks and is defined by 

VV: mi = m2 = xm^ei, m^''^ = m^ei, (12) 

which leads to the expressions 

R^{x) = l--J-i/rcf[(2x-l)lnx + 2(x-l)lnyrcf] 

iVf 

-l/ref(a; - 1) [(2a8 - "5) + ^] (13) 

Rnx) = 1 - ^2/rcf [{x + 1) In (i (x + 1)) + (x - 1) In y,,f 

+yrcf{x - l)la5, (14) 

where A^f denotes the number of dynamical quark flavours. The case labelled "VSl", 
based on non-degenerate valence quarks, is defined by 

VSl: nil = xnirei, m2 = rrf^^ = mref- (15) 



According to Table 1 of ref. |0 the expressions for the ratios are then given by 

^M^^(^) = 1 - ^T^Z/rcf [a;lnx + (a;- l)ln?/i.ef] -l/ref(a;- I)|(2a8-a5) (16) 

iVf 



Rr\x) = 1-^y.ef 



[x + 1) In (|(x + 1)) + {x-l) Iny^ef + ^\nx 



+y^^,{x-l)l[a, + ^). (17) 

Our simulations were performed for Nf = 2 flavours of dynamical, 0(a) improved Wilson 
fermions. The value of the bare coupling was set to (3 = Q/Qq = 5.2. The improvement 
coefficient Csw, which multiplies the Sheikholeslami-Wohlert term in the fermionic action, 
was taken from the interpolating formula of ref. |T^ . Here we considered a single value 
of the sea quark mass, corresponding to a hopping parameter k,^^'^ = 0.1355. For this 
choice of parameters we generated 208 dynamical gauge configurations on a lattice of size 
16^ ■ 32. For further details we refer to 0-18]. Here we only mention that the hadronic 



radius tq defined through the force between static sources ^9j has been determined as 
To/a = 5.041(40) [|I8|. For tq = 0.5 fm this implies that the lattice spacing in physical 
units is a = 0.099(1) fm. 

We have computed quark propagators for valence quarks with hopping parameters k™' = 
0.1340, 0.1345, 0.1350, 0.1355 and 0.1358. In addition to local operators, we have also 
calculated propagators for fuzzed sinks and/or sources, using the procedure of [^. In 
the pseudoscalar channel we employed both the pseudoscalar density and the temporal 
component of the axial current as interpolating operators. These two types were used 
together with the different combinations of fuzzed and local propagators to construct a 
4x4 matrix correlator for pseudoscalar mesons. By performing factorising fits using an 
ansatz that incorporates the ground state and the first excitation, we were able to extract 
the pseudoscalar mass mps, as well as the matrix elements of the axial current and the 
pseudoscalar density, i.e. 

Ca = (0|Ao|PS), Cp = (0|P|PS). (18) 

In order to be consistent with 0(a) improvement the amplitudes Ca and Cp must be related 
to the matrix elements of the improved currents and densities at non-zero quark mass. 
Using the definitions of [^ it is then easy to see that the (unrenormalised) pseudoscalar 
decay constant is given by 



/Ca ^^ c 



Fps = (1 + hxam^) I ^^^ + ca^^ sinh(amps) ) , (19) 

( mps mps J 

whereas Gps is given by 

Gps = (1 + 6pamq)Cp- (20) 

The improvement coefficients 6a , ^p and ca were computed in one-loop perturbation 



theory p2|,|23| in the bare coupling, since non-perturbative estimates are not available at 



present. For the improvement coefficient ca this procedure yields 

ca = -0.0087. (21) 

For degenerate valence quarks the quark mass m^ is given by 

where we have inserted Kcrit = 0.13693 for the critical hopping parameter, as estimated 
ISf . Finally we note that the current quark mass m defined through the PCAC relation 



m 



in the 0(a) improved theory is obtained from ^a, Cp ^^^ "^ps via 

am = —— sinh(a?Tips) H — CAsinh^(amps). (23) 

2Cp 2 

All our statistical errors were obtained using a bootstrap procedure 



4 Results 

Our results for pseudoscalar masses, matrix elements and current quark masses are listed 
in Table ^ Compared with [|I^ the numbers for the pseudoscalar masses reported here 
may differ by up to one standard deviation, as a result of using a larger matrix correlator 
in the fitting procedure. 

Following eq. (0) we have defined the mass rriref at the reference point for k™' = /t™' = 
^sea _ 0.1355, which corresponds to (^o"^Ps)^lm=m ~ 2.092. Using the leading-order 
relations m|g = 23^171 and Fq = F^^ = 93.3 MeV this implies 

l/ref = 0.2370. (24) 

It is instructive to compare these values with those of the previous quenched study 0, 
where ?/ref = 0.3398, (romps)^|m,=mrof = 3, and mref ~ rris. Hence, the value of m-^ei 
employed in this paper is smaller by 30%, such that rriref ~ 0.7r7is. A smaller value 
of TTiref IS clcarly desirable for our purpose, since the predictions of ChPT are expected 
to hold more firmly for smaller masses. Nevertheless, our sea quarks are still relatively 
heavy, as signified by the ratio {jn-ps / m\)\m=m^<^^ ~ 0.58 [l^, which is to be compared to 
m.„/mp = 0.169. 

In the spirit of partially quenched QCD we have considered valence quarks that are 
lighter than the sea quarks. We were thus able to extend the quark mass range down 
to about ms/2, which is quite a bit below the smallest mass reached in 0. This may 
be an indication that the inclusion of dynamical quarks alleviates the problem of excep- 
tional configurations, which precludes attempts to work at very small quark masses in 
the quenched approximation. However, our attempts to push to valence quark masses be- 
low ms/2 have proved unsuccessful, due to the appearance of exceptional configurations 
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^val 
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amps 


aCA/mps 


a'Ci 
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am 




0.1358 


0.1358 


0.2301 


+ 52 
-51 


0.0995 


+ 24 
-30 


0.1782 


+70 
-57 


0.0147 


+ 4 
-5 


0.1355 


0.1355 


0.2869 


+40 
-41 


0.1068 


+ 22 
-25 


0.1893 


+57 
-52 


0.0232 


+ 4 
-5 


0.1350 


0.1350 


0.3585 


+ 26 
-28 


0.1166 


+ 20 

-24 


0.1046 


+47 
-49 


0.0368 


+ 4 
-5 


0.1345 


0.1345 


0.4192 


+ 20 
-23 


0.1246 


+ 19 
-24 


0.2187 


+40 
-45 


0.0507 


+ 5 
-5 


0.1340 


0.1340 


0.4739 


+ 17 
-20 


0.1312 


+ 17 
-25 


0.2315 


+43 
-44 


0.0650 


+ 6 
-5 


0.1358 


0.1355 


0.2607 


+45 
-46 


0.1033 


+ 24 
-26 


0.1841 


+60 
-55 


0.0190 


+ 4 
-4 


0.1350 


0.1355 


0.3249 


+30 
-37 


0.1118 


+ 20 
-25 


0.1970 


+50 
-52 


0.0300 


+ 4 
-5 


0.1345 


0.1355 


0.3591 


+ 27 
-32 


0.1158 


+ 20 

-24 


0.2036 


+50 
-46 


0.0369 


+ 4 
-5 


0.1340 


0.1355 


0.3907 


+ 25 
-28 


0.1190 


+ 21 
-25 


0.2094 


+44 
-45 


0.0438 


+ 4 
-6 



Table 1: Results for pseudoscalar masses, matrix elements and current quark masses at /? = 5.2, 

K^ea ^ 0.1355. 

for TT^vayj^sca < q^ Moic pieclsely, we observed large statistical fluctuations in hadron 
correlators computed from quark propagators at k™^ = 0.1360, despite the fact that the 
inversion algorithm converged. 



The results in Table |] can now be used to compute Ru and Rp through eqs. (|3), (|19D 



and (0). In order to map out their quark mass dependence in some detail, we have 
performed local interpolations of the results for Fps and -Fps/Gps to 20 different values of 
the dimensionless mass parameter x in the range 0.7 < x < 2.6, separated by increments 
of 0.1. 

Since we only have one /5- value we cannot extrapolate Ru, Rp to the continuum limit for 
fixed X. Unlike ref. 0, where such extrapolations could be performed, we must compare 
our data to ChPT at non-zero lattice spacing. Our estimates for the low-energy constants 
are therefore subject to an unknown discretisation error. It is reasonable to assume, 
however, that cutoff effects in Rm and Rp are fairly small, owing to cancellations of 
lattice artefacts of similar size between numerator and denominator. We will return to 
this issue below, when we discuss our final estimates for the low-energy constants. 

In order to extract the low-energy constants we have restricted the x-interval to 0.7 < x < 
1.1, thereby seeking to maximise the overlap with the domain of applicability of ChPT, 
whilst maintaining a large enough interval to check the stability of our results. Estimates 
for the low-energy constants were obtained by fitting the data for Rm and Rp to the 
corresponding expressions for the "VV" and "VSl" cases listed in Section ^. Since Rp is 
linear in a^ one can obtain this low-energy constant also from simple algebraic expressions 
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(b): Rm 
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Figure 1: (a): The ratios Rp for the cases labelled "VV" and "VSl" compared with the fit to 
the expressions in ChPT; (b): the same for the ratios Rm- Dotted lines indicate the variation 
due to the statistical uncertainty in the low-energy constants. The straight double lines indicate 
the quark mass range explored in the earlier quenched study 0]. 



involving the difference Rf{xi) — Rf{x2) for two distinct arguments, xi and X2- A similar 
relation can be used to compute {2as — as) from Rm{xi) — Rm{x2)- We have checked that 
both methods give consistent results and obtain 



W: 

VSl: 



(2) 



a 



(2) 



1.20 



1.22 



+11 

-16' 



-11 
-16' 



(2a 



(2) 



a 



(2)^ 



(2ar-«r: 



0.36 
0.36 



-10 
-10' 

-10 
-12' 



(25) 
(26) 



where the errors are purely statistical. From here on we also indicate the number of 
dynamical quark flavours as a superscript, to distinguish these estimates from the corre- 
sponding ones in the quenched and three-flavour cases. Our results can now be inserted 
back into the expressions for Rm and Rp. The resulting curves are plotted together with 
the data in Fig. |l|. 

It is striking that the data for Rp in the VV case are described remarkably well over 

(2) 

the whole mass range, despite the fact that al has only been determined for a; < 1.1. 
The qualitative behaviour of Rp - which features a slight curvature - is thus rather 
well modelled by eqs. (|1^ and (p!7D, which include a linear term as well as chiral loga- 
rithms. By contrast, there are no logarithmic contributions to Rp in quenched ChPT, 
and the expected purely linear behaviour has indeed been observed in the data 0. Of 
course, higher-order terms in the quark mass could in principle produce a curvature, and 
therefore these observations do not provide unambiguous evidence for chiral logarithms. 
Nevertheless, it is remarkable that the clear distinction between the expressions for Rp 
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in partially quenched and quenched ChPT (i.e. the presence, respectively absence of chi- 
ral logarithms) is accompanied by corresponding qualitative differences in the numerical 
data. 

Unlike Rp the ratio Ru is not described well for larger masses, which may signal a break- 
down of the chiral expansion for this quantity for masses not much larger than rris. It is 
therefore conceivable that higher orders in ChPT affect the extraction of (2a8 — a^ ). 
However, without access to quark masses that are substantially lower than our simulated 
ones, it is not easy to quantify reliably any uncertainty due to neglecting higher orders. 

One way to examine the influence of higher orders is to extract the low-energy constants 
from a mass interval of fixed length, (a;max ~ ^^min) = 0.4, which is then shifted inside 
an extended range of 0.7 < x < 1.5. The spread of results so obtained then serves as 
an estimate of the systematic errors incurred by neglecting higher orders. We note that 

(2) 

X = 1.5 corresponds to a quark mass slightly larger than mg. For 0:5 such a procedure 
yields only a small variation of ±0.05. This is not surprising, since Rp is modelled very well 
over the entire mass range. By contrast, the spread of results obtained for (2ag — al ) 
is as large as ±0.15. Whether or not these numbers represent realistic estimates of the 
actual uncertainty cannot be decided at this stage. In order to be more conservative 
we have decided to quote a systematic error of ±0.2 for all low-energy constants. We 
note that this level of uncertainty due to neglecting higher orders was also quoted in the 
quenched case |§], where quark masses were shghtly larger. 

Since we do not have enough data to extrapolate Rp and -Rm to the continuum limit 
we also have to estimate a systematic error due to cutoff effects. As explained above, 
however, we expect such effects to be small. In order to get an idea of the typical size of 
discretisation errors we have looked again at quenched data obtained at (3 = 5.93 |I6HT8[ 



and 6.0 [^, for which the lattice spacing in physical units is roughly the same as in 
our dynamical simulations (a ~ 0.1 fm). For both (3 = 5.93 and 6.0 the results for 
Ru and Rp are mostly consistent within errors with the corresponding values in the 
continuum limit (see, for instance. Fig. 1 in [^]). Furthermore, low-energy constants 
extracted for a ^ 0.1 fm differ from the results in the continuum limit by less than 
one standard deviation. Although these findings cannot be taken over literally to the 
dynamical case without direct verification, they nevertheless indicate that lattice artefacts 
are small enough such that does not have to expect large distortions in our estimates for 
the low-energy constants. In order to take account of these observations we have decided 
to quote an additional systematic error due to lattice artefacts, which is as large as the 
statistical error. 

Since non-perturbative estimates for the improvement coefficients b\, bp and ca are not 
available for A^f = 2, one may be worried that there are large uncancelled lattice artefacts 
of order a in our data. We have addressed this issue by studying the influence of different 
choices for improvement coefficients on our results. To this end we have repeated the com- 



plete analysis using non-perturbative values for ca and the combination 6a — ^p obtained 
in the quenched approximation ||25|,^ at a similar value of the lattice spacing, a ~ 0.1 fm. 
We found that the resulting variation in the estimates for a^ and {2ag — 0:5 ) is typi- 
cally a factor 10 smaller than the statistical error. Thus we conclude that the influence 
of improvement coefficients on our results is very weak indeed. 

As a final comment we point out that we have not taken finite volume effects into account 
in our error estimates, because different lattice sizes were not considered in our study 



(unlike in earlier simulations ||2^). However, since Lmps = 4.59 at the reference point 
and Lrrips = 3.68 at the lightest valence quark mass, one may not be totally convinced 
that such effects may be entirely neglected. We stress, though, the the definition of Rm 
and Rp implies that only the relative finite-size effects between hadronic quantities is 
relevant. Thus, as long as the mass parameter x does not differ too much from unity, 
one can reasonably expect that finite-volume effects largely cancel in the ratios used to 
determine the low-energy constants. Finite-size effects in pseudoscalar masses and decay 



constants have also been studied in ChPT 28-311. These calculations indicate that the 



typical relative finite- volume effect in mps and Fps between our reference point and the 
smallest quark mass is less than 1%. 



5 Discussion and outlook 

After combining the different systematic errors in quadrature we obtain as our final results 
in two- flavour QCD: 

af = 1.22+°:}^ (Stat) +0:1 (syst) (27) 

(2af^-af^) = 0.36 ± 0.10 (stat) ±0.22 (syst) (28) 

af^ = 0.79+°:°^ (stat) ±0.21 (syst), (29) 

where eqs. (|27D and ( pH]) have been combined to produce the result for ag . 

We can now investigate the dependence of the low-energy constants on the number of 
dynamical quark flavours. In the quenched approximation (i.e. for A^f = 0) it was 
found thatQ 

aP = 0.99 ± 0.06 (stat) ± 0.2 (syst) (30) 

aP = 0.67 ±0.04 (stat) ±0.2 (syst). (31) 

A comparison with eq. (P7| ) and (^) then shows that al and Og are larger than their 
quenched counterparts by 23% and 18% respectively. We can thus conclude that the Af- 
dependence of the low-energy constants is fairly weak: variations between the quenched 
and two-flavour theories are about as large as the error due to neglecting higher orders. 

^When extracting the result for a^ it was assumed that the coefScients multiplying quenched chiral 
logarithms were set to S = 0.12, a* = 0.0. 
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Figure 2: Data for Ryi (VV case) and the curve which results if A'^f = 3 in the determination 
of (2ag — a^ ) (solid line). Dotted lines represent the statistical uncertainty. The dashed 

(3) (3- 



(3) (3) 

curve corresponds to (2ag — 05 ) = —2.3, a value which is consistent with the hypothesis of a 



massless up-quark. 



Although there is a priori no reason why the weak A/'f-dependence should extend to the 
physical three-flavour case, it is still instructive to compare eqs. (|27| ) and ( pOf ) with phe- 
nomenological values of the low-energy constants. It then becomes obvious from eqs. (H) 



and (0) that our results for ai, and «§ are compatible with the standard estimates 
found in the literature. By contrast, our numerical data for Rp and Rm suggest that a 
large negative value for ag, which is required for the scenario of m^ = (see eq. (||)), 
is practically ruled out. Thus, provided that the quark mass behaviour in the physical 
three-flavour case is not fundamentally different, the possibility of a massless up-quark is 
strongly disfavoured. 

By how much does one expect the mass dependence of Rm and Rp to differ between Nf = 2 
and 3 ? Ultimately this must be answered by a direct simulation of the three-flavour case. 
For the time being we have to be content with the following gedanken simulation. Suppose 
that we had analysed our N{ = 2 data under the erroneous assumption that they had been 
obtained in the physical three-flavour case. We would then have set A^f 



and ( [T7| ) to extract a 



(3) 



• • (3) 

givmg al 



0.98 



+ 11 +23 
-16 -26' 



3 in eqs. (|1^ 
This value can be inserted into the 



expression for Fk/-FV in the physical theory [TTI], to yield 



— = 1.247" 



-0.013 (statj _o.o2i (systj. 



(32) 



which is in fair agreement with the experimental result Fk/Ftt = 1.22 ± 0.01. This shows 
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that the experimental value can only be reproduced if the quark mass dependence of Rp 
in the physical case is not much different from that encountered in our simulations for 
Nf = 2. In other words, it is reasonable to assume that the mass dependence of Rp is 
only weakly distorted by neglecting the dynamical quark effects due to a third flavour. 

Similarly we can apply the expressions for Rm for A^f = 3 to our data, which gives 

-10 +22 
-12 -23' 

mass correction Am is then obtained as 



(2o;g — 0:5 ) = 0.20 ^-,„ ^„„. The corresponding curve is shown in Fig. El. The first-order 



,^ = -0.041°:°^ (stat) ± 0.11 (syst), (33) 



which is consistent with Leutwyler's estimate (see eq. (^). We emphasise that this does 
not represent a reliable result for Am derived from first principles. Nevertheless, the above 
discussion shows that there are examples which support the idea that the gross features 
of the mass dependence do not differ substantially in the two- and three-flavour cases. 
On the basis of this assumption one may conclude that the correction factor Am is indeed 
small, ruling out the scenario of a massless up-quark. As a further illustration we have 
included in Fig. |^ the curve which one would expect if ^m^^ = 0, by taking the central 
values for a^ and ag from eqs. (0) and (|^) . 

The first priority for future work is undoubtedly the application of the method to simula- 
tions employing A^f = 3 flavours of dynamical quarks, and the extension of the quark mass 
range towards the chiral regime. While efficient simulations of QCD with odd A^f and light 
dynamical quarks represent an algorithmic challenge, some efforts in this direction have 
already been made [0. It would also be interesting to extend applications to the case of 



flavour singlets, which allow a determination of 07 P], i.e. another low-energy constant 
afflicted with the KM ambiguity. Methods to improve the notoriously bad signal/noise 
ratio in flavour-singlet correlators have been developed [^^, so that there are good 
prospects for a successful implementation. 
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